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Abstract 

A lattice model of critical dense polymers is solved exactly for arbitrary system size 
on the torus. More generally, an infinite family of lattice loop models is studied on 
the torus and related to the corresponding Fortuin-Kasteleyn random cluster models. 
Starting with a cylinder, the commuting periodic single-row transfer matrices are built 
from the periodic Temperley-Lieb algebra extended by the shift operators f^ 1 . In 
this enlarged algebra, the non-contractible loop fugacity is a and the contractible loop 
fugacity is (5. The torus is formed by gluing the top and bottom of the cylinder. This 
gives rise to a variety of non-contractible loops winding around the torus. Because 
of their nonlocal nature, the standard matrix trace does not produce the proper ge- 
ometric torus. Instead, we introduce a modified matrix trace for this purpose. This 
is achieved by using a representation of the enlarged periodic Temperley-Lieb algebra 
with a parameter v that keeps track of the winding of defects on the cylinder. The 
transfer matrix representatives and their eigenvalues thus depend on v. The modified 
trace is constructed as a linear functional on planar connectivity diagrams in terms of 
matrix traces Tr^ (with a fixed number of defects d) and Chebyshev polynomials of the 
first kind. For critical dense polymers, where (3 = 0, the transfer matrix eigenvalues 
are obtained by solving a functional equation in the form of an inversion identity. The 
solution depends on d and is subject to selection rules which we prove. Simplifications 
occur if all non-contractible loop fugacities are set to a = 2 in which case the traces are 
evaluated at v = 1. In the continuum scaling limit, the corresponding conformal torus 
partition function obtained from finite-size corrections agrees with the known modular 
invariant partition function of symplectic fermions. 
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1 Introduction 



Solvable critical dense polymers CAi(l, 2) is the first member of the Yang-Baxter integrable 
family £A4(p,p') of logarithmic minimal models [1]. Algebraically, the logarithmic minimal 
models are described by the planar Temperley-Lieb (TL) algebra [2, 3] with loop fugacity 
8 = 2 cos A and crossing parameter A = iH_z£lZL. p or solvable critical dense polymers, A = f 
so the loop fugacity vanishes, 6 = 0. This implies that closed contractible loops are not 
allowed so the loop segments form long polymer segments. An elementary face of the square 
lattice is assigned a statistical face weight according to the configuration of the face. The 
two possible configurations with their associated weights are combined into a single face 
operator as 



u 



sin (A 



u 



+ sinu 



A- 



where u is the spectral parameter and the lower-left corner has been marked to fix the 
orientation of the square. By the crossing symmetry, rotating the face by 90 degrees changes 
the spectral parameter from u to A — u. The polymer segments begin and end at nodes at 
the midpoints of the edges of the face. 

Because of the nonlocal degrees of freedom, in the form of long polymer segments, the 
topology of the lattice has profound effects on the properties of dense polymers. Critical 
dense polymers has been solved exactly for arbitrary finite sizes on the strip [4-6] and the 
cylinder [7]. Among other results, these studies have firmly established the central charge 



-2 and an infinitely extended Kac table of conformal dimensions A r 



r, s 



1,2,3. 



given by the usual Kac formula. The integrals of motion and Baxter Q-operators have also 
been studied [8]. Recently, generalized order parameters were calculated [9] for an exactly 
solvable <p\$ off-critical perturbation of dense polymers. For (r, s) satisfying (2r — s) 2 < 8, the 
conformal dimensions A r s = 8 TjS precisely agree with the critical exponents 6 r>s associated 
with these generalized order parameters. 

Perhaps, the only other Yang-Baxter integrable model that has been studied in such 
detail is the Ising or free fermion model. The glaring anomaly is that the lattice model 
of critical dense polymers has not yet been solved exactly on the torus. The fundamental 
obstacle in applying transfer matrix techniques is that, because of the nonlocal degrees of 
freedom and winding, the naive matrix trace of a cylinder transfer matrix does not produce 
the proper geometric torus [7]. In this paper, we use link representations that depend on a 
winding parameter v and construct a modified trace, initially as a linear functional on planar 
connectivity diagrams, and ultimately in terms of matrix traces Tr^ (with a fixed number of 
defects d). These representations, which we refer to as twist representations and denote by 
Ud, were first studied in [10, 11] and more recently in [12]. 

Modified matrix traces have appeared before as a way to relate Potts models and loop 
transfer matrices [13, 14]. In [13], the partition function of the Potts model on the strip is 
expressed in terms of representations labeled by a number I of bridges, isomorphic to link 
state representations with I defects. The coefficients of the corresponding modified matrix 
trace are given by Chebyshev polynomials of the second kind. On the torus, this construction 
fails, but an alternative construction was subsequently obtained by the same authors [14] 
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using a larger class of representations involving the cyclic group Cg. The coefficients of the 
corresponding modified matrix trace are given by Chebyshev polynomials of the first kind. 

Here we use a different construction based on the twist representations u^, and to the 
best of our knowledge, the ensuing modified matrix trace is new. Furthermore, the fugacities 
of non-contractible loops considered in [13, 14] are very specific, while the ones we consider 
are general. In principle, for critical dense polymers, this allows us to calculate the torus 
partition function exactly for arbitrary finite system sizes and arbitrary non-contractible 
loop fugacities by solving an inversion identity in the form of a functional equation for the 
transfer matrix eigenvalues in sectors with d defects. 

In general, the double-row transfer matrices on the strip [1, 15] and the periodic single- 
row transfer matrix on the cylinder [7, 16] exhibit nontrivial Jordan cells. In this paper, 
we are interested in partition functions so we only study the eigenvalues, conformal spectra 
and associated characters, and do not consider the structure and indecomposability of the 
associated Virasoro representations. 

Using the modified trace, and calculating finite-size corrections, allows us to obtain the 
torus partition function in the continuum scaling limit. From general principles of conformal 
invariance [17], this is expected to be modular invariant. Indeed, we find that the torus 
partition function is Coulombic [18], modular invariant and in agreement with [19]. This 
modular invariant partition function coincides with that of the triplet model [20] and sym- 
plectic fermions on a Z 2 orbifold [21]. Both of these models are described by a logarithmic 
Conformal Field Theory (CFT) [22] with central charge c = -2 [23-29]. 

A striking feature of our derivation of the modular invariant is that, contrary to most 
results on loop models and the associated conformal field theories, it is entirely rigorous. In 
particular, we present a mathematical proof of the selection rules determining the eigenvalues 
(and their degeneracies) of the transfer matrix. This proof is based on a relation between the 
loop model and an XXZ Hamiltonian, an idea previously used in [5] to construct a similar 
proof of the selection rules applicable to the loop model defined on the strip. A proof of the 
selection rules conjectured in [7] is provided in the process. 

The layout of this paper is as follows. In Section 2, we define a general loop model on 
the torus and explain its relation to the Fortuin-Kasteleyn (FK) random cluster model [30]. 
We also introduce the enlarged periodic Temperley-Lieb (TL) algebra. Lastly, we introduce 
link states, the twist representation on these link states, which uses the winding parameter 
v to keep track of windings, and our new modified trace which we use to close the cylinder 
to the geometric torus. In Section 3, we focus on critical dense polymers. We use the twist 
representation in sectors with d defects and the known transfer matrix inversion identity [7] to 
write down functional equations involving v for the eigenvalues. These equations are solved 
for the eigenvalues subject to certain selection rules which we prove in Appendix A. Our 
modified trace is then applied to build the finitized torus partition function. In the continuum 
scaling limit, this becomes a simple Coulombic partition function and so is modular invariant. 
Loop models and critical dense polymers on helical tori [31-33] are considered in Appendix B. 
Finally, Section 4 contains some concluding remarks. 
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2 Loops, FK clusters and the periodic TL algebra 



2.1 Loop model 

Configurations of the loop model of our interest are drawn on a square lattice of M x 
N tiles in the plane. Every tile is decorated with one of the two diagrams 



and 



as indicated in Figure 1. The boundary conditions are taken to be periodic in both the 
vertical and horizontal directions. The resulting diagram is a graph of non-intersecting 
closed curves (loops) onanMxAf torus (of trivial helicity, see Appendix B). Because of the 
toroidal geometry, the loops can have different homologies. Contractible loops, which can 
be continuously deformed to a point, are said to be homotopic to a point, or of homotopy 
{0}. Loops can also wind around the torus in a nontrivial manner; a times in the horizontal 
direction and b times in the vertical direction. The constraint that the loops do not intersect 
imposes that the greatest common divisor of a and b, a A b, is 1 (with aA0 = 0Aa = a). 

As a curve moves upward to wind around the torus vertically (b times), it can wind around 
the horizontal direction (a times) by either going to the right or the left. To establish the 
difference in homotopy between these two cases, we fix the convention that a is positive if 
the curve winds toward the right, and negative if it winds toward the left. The ensuing loop 
will be said to have homotopy {a, b}, with o 6 Z, ft 6 N and a A b = 1. For example, the 
configuration in Figure 1 has two loops of homotopy {1, 1}. 



M 




N 



Figure 1: A configuration of the loop model with weight Wl(ct) = f3 A a\ l sin 8 wsin 16 (A — u). 



In a given configuration, there can be more than one loop that winds nontrivially around 
the torus, but because they are non-intersecting, these non-contractible loops must all have 
the same homotopy. A configuration is therefore classified by the homotopy of its loops 
that wind around the torus. If a configuration has no loop winding around the torus, it has 
homotopy {0}. 

To every configuration a, we associate a weight 

WM = r(\{a^y i { ®p n 2 m (2.1) 

aAb=l 



where n(E2) and n(s) are the number of times the tiles YA and appear in the configura- 
tion a, Uq is the number of loops of homotopy {0} and n a ^ the number of loops of homotopy 
{a, b}, while riaA6=i * s a P r °duct over a G Z and b G N with the constraint a A b = 1. Note 
that at most one n a ^ is nonzero, so the product will contain at most one term ct a a b b . If no 
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loop winds around the cylinder, then the product is replaced by 1. We refer to (3 and a a ^ 
as contractible and non-contractible loop fugacities, respectively, while pi and P2 are the 
weights of the face configurations ^ and respectively. 

The partition function is defined as the sum of the weights of the 2 MN configurations, 



5> 



l{o- 



(2.2) 



From the discussion above, the sum in can be split between the different possible ho- 
motopies of the configurations, and we write Z^ih) for the partition function obtained from 
restricting the sum over a to configurations with homotopy h. We thus have 



Z L = Z L ({0})+ E Z d{*M), 



(2.3) 



a A 6=1 



where ^ aAfe=1 is a sum over a £ Z and &GN with the constraint aAb= 1. By construction, 
for M and N finite, Zi({a, b}) is zero if \a\ > M or b > N. 

The partition function depends on many free parameters: (3, a a fi, p\ and p 2 . Since 
ft (J2) + = MiV, up to a cx-independent constant, Wl(ct) depends only on the ratio of 
Pi and p2- In the following, we parameterize the weights Pi,P2 and the contractible loop 
fugacity (3 as 

p 1 = sm(\ — u), p 2 = sinw, /3 = 2cosA, (2.4) 
where A and u are the crossing and spectral parameters, respectively. 



2.2 Relating the loop and FK cluster models 

The Fortuin-Kasteleyn (FK) random cluster model [30] can also be defined on the toroidal 
geometry of the loop model above. We label the lower-left corners of the tiles by their 
positions (x,y) on Z 2 , with (0,0) in the lower-left corner of the lattice. A distinguished 
sublattice is added to the lattice: Sites, identified by solid dots in Figure 2, appear in the 
positions (x, y) if x + y is odd. The integers M and N are chosen even to ensure periodicity 
in both directions. As before, a configuration consists of a choice of either ^ or ^ for each 



M 




N 

Figure 2: An FK configuration with weight Wfk{&) = Q^W\^v\v2 ■ 



of the tiles. Dashed lines are then added between neighbouring sites of the distinguished 
sublattice if the configuration allows one to draw a bond without intersecting the loop 
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segments. Blue lines ( #**) are for bonds connecting upper-right and lower-left corners of the 
faces (type 1 bonds), while red lines (\ ) are for bonds connecting upper-left and lower- right 
corners of the faces (type 2 bonds). For a given configuration, we denote by N Bl and N B2 
the numbers of type 1 and type 2 bonds. 

The curves formed by concatenating the loop segments of the tiles separate the set of 
sites into clusters consisting of sites connected by bonds. The definition of a cluster site is 
therefore intimately related to the notion of the interior of a closed loop, here defined to 
be occupied by sites. For M or N odd, the periodicity of the lattice renders this definition 
ill-defined, which, as already announced, prompts us to only consider FK clusters on lattices 
for which both M and N are even. 

As for loops, homology properties are well-defined for these clusters, some of which are 
contractible and said to have cluster homotopy {0}. The contour of such a cluster is a 
union of loops with homotopy {0}. Other clusters have homotopy {a, b} and wind around 
the cylinder a times in the horizontal direction and b times in the vertical direction (with 
a A b = 1). The boundary of such a cluster contains two loops of homotopy {a, b}. Finally, 
unlike loops, a cluster of sites can wind around the torus in another nontrivial manner: by 
wrapping it in both directions. Such a cluster appears in the first configuration in Figure 3. 
Such clusters are said to have cross-topology, or homotopy Z 2 . Only a single cross-topology 
cluster can appear in a given configuration, and it cannot coexist with any cluster with 
homotopy {a, b}. The boundary of a cluster with cross-topology consists of loops with 
homotopy {0}. By construction, FK configurations are characterized by the homotopy of 
their clusters that wind the torus nontrivially. 

The FK weight of a configuration is 

W FK (a) = Q»° ( J] <r>Nf Sl ^ 2 (2-5) 

aAb=l 

where Q, w a ,b, w + , V\ and v% are all free parameters, while N , N atb , N + are the respective 
numbers of clusters homotopic to a point, of homotopy {a, b}, and with cross-topology. Only 
N + or one N a>b can be nonzero, and N + G {0, 1}. The partition function Z FK is the sum of 
Wfk(o~) over all possible configurations and can be separated into contributions coming from 
each homotopy group. We thus write cr(h) for the set of configurations a with homotopy h, 
and Yla(h) f° r ^ ne sum over configurations with homotopy h. 

Our next objective is to show how to calculate Zfk at the critical point from Zl, the 
partition function of the loop model. To proceed, we note that there exists an Euler relation 
between the cluster numbers N and N a)b , the total number of bonds N B = N Bl + N B2 , the 
number of sites N s = NM/2 and the number of loops in the diagram, #(cr) = n + n a y. 

N + N a , b = \ + N s - N B ) . (2.6) 

This can be verified by first checking that it holds for the configuration do with zero bonds, 
where N = #(cr ) = NM/2 and N ajb is zero. By adding bonds, one either reduces the 
number of clusters, adds cycles in clusters, creates {a, b} clusters or creates a Z 2 cluster, 
and equation (2.6) remains satisfied in all cases. We note that the number of cross-topology 
clusters, N + , does not appear in the relation (2.6). 
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With these observations, we can write the FK partition function as 



(\ N a b 

m ;»! I n n .,. ,7, J ^ ' <j(Z 2 ) 

E + E E )<>»<*<**> f^"" f^"" 

cr({0}) aAb=l <t({o,6}) 



Q J \VQJ \VQ 



a(Z 2 ) 



, Nb, / \ Nfl, 



+ « + (-^j (, 7) 



= Q |Ar *(E + E E +-+ E)^ |no K^) |n<a (-^) 

CT ({0}) aAfc=l er({o,6}) a(Z 2 ) WW 

where the equation (2.6) and the relation N a ^ = §?V& were used in the second equality. 
Furthermore, our notation uses implicitly that, for configurations of homotopy {0}, N a ^ = 
for all a A b = 1 . 

Two more ingredients are required to proceed further. First, we set the parameters v\ 
and t>2 to critical values [34] by imposing the relation 

v x v 2 = Q. (2.8) 

These values are parameterized by the crossing parameter A and the spectral parameter u 
through the relations 

vi VQ sin(A-u) r- 

\JQ = 2 cos A. (2.9) 



y/Q v 2 sinii 

The other key element is a bijection that exists between configurations with cross- 
topology and configurations with homotopy {0}, obtained by shifting every tile one position 
to the right. An example is given in Figure 3. (Equivalently, one could interchange the 
distinguished and non-distinguished sublattices.) The number n Q is the same in the two 
configurations, whereas the numbers of bonds typically differ. The bonds in the first con- 
figuration are all absent in the second one, and vice versa. For i = 1,2, let N Bi and N' B . 
denote the numbers of type % bonds in the configuration with cross-topology and homotopy 
{0}, respectively. It then follows that 

N Bl + N' B2 = NM/2 = N' Bl + N B2 (2.10) 

and hence 

N Bl - N B2 = N' Bl - N' B2 . (2.11) 
From the above discussion, we now find 



J FK 



9T1({1 + W+ )Z({0})+ E ^({ fl ' & }))> (sin M sin(A- M ))^ M (2.12) 



a A 6=1 
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Figure 3: An illustration of the bijection between a cross-topology configuration and a 
configuration with homotopy {0}. The second configuration of tiles is a shift towards the 
right of the first. 



where 

Z({0}) = Q^ no (smu^ NM+NB ^ N ^ (sin(A - u ))* nm+Nb i- Nb * , (2.13) 
Z({a,b}) = Q^(«^)^ k (sintO*™^ (2.14) 

*({a,b}) 

To express these in terms of the partition function of the loop model Zl, we note that the 
powers of sinw and sin(A — u) are given by 

n(H) = \NM + N B2 - N Bl , n($) = \NM + iV Bl - N B , 2 . (2.15) 

By setting /3 = \/Q and a a ,b = y/u>a,b, we find Zi({a,b}) = Z({a,b}). The sum over 
configurations in Zl({0}) of the loop model contains all configurations with only loops that 
are homotopic to a point, and is therefore ^ CT ({ }) +J2a(z 2 )- By setting (3 = \[Q and using 
the previous bijection, we find Z L ({0}) = 2Z({0}). Finally, setting all non-contractible loop 
fugacities a a ^ to yields Z L \ a b=0 = Z L ({0}). In conclusion, we have found that the FK 
partition function can be written as 



Q-2 N ° w+-l 



Z, 



«a,6=0 



(2.16) 



2.3 FK model with Q = 

Even though the FK and loop models both are well-defined in the /3 — > limit, the passage 
from the first to the second is awkward for /3 = 0: The factor i n equation (2.16) goes 

to zero, the critical condition requires that either v\ or t>2 (or both v\ and V2) be zero, and 
the parameterization (2.9) of v\ and v 2 in terms of the variable u is ill-defined. Here we show 
how to treat the case y/Q — (3 — properly. 

In the FK model, setting Q = is equivalent to giving a weight zero to configurations 
with clusters of homotopy {0}. We implement this by writing 

W FK {a) = 5 Nofi ( J] <r)^Nf fll ^ fl2 (2-17) 

aA6=l 
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and accordingly 

Z ™=zZ E S Not0 w^ B ^+w + J2S No ,oV^v^. (2.18) 

a/\b=l a({a,b}) 



To link the FK model with Q = to the corresponding loop model, we introduce a restricted 
FK model by imposing the further constraint on the configurations that 

the clusters contain no contractible cycles of bonds. 

This step is crucial. In the original FK model, loops with homotopy {0} occur in three 
scenarios: They can act as a boundary for cross-topology clusters or for {0} clusters, or 
appear if a non-contractible cluster contains a contractible cycle of bonds. With the added 
no-cycle constraint, the last two scenarios are prohibited. One can therefore replace £tv ,o by 
a similar constraint on no, the number of loops of homotopy {0}: no = if a G a({a,b}) 
but no = 1 if o G o~(Z 2 ). The corresponding restricted partition function is denoted by Z' FK 
and, from the bijection between {0} and Z 2 configurations, is given by 



ryl r N Bi N B ST-^ X l n a,f> N Bi N B 2 

Z FK =W + 2^<W U 1 V 2 + }^ 1^ 5 »ofl W lb V l V 1 

a/\b=l a({a,b}) 



= £j (t E + E E ) «f . (a w) 

<r({0}) aAfe=l cr({a,6}) 

where we again use implicitly that n a ^ = for all configurations of homotopy {0}. The limit 
P in (2.19) is well-defined. 

To relate this restricted FK model to the loop model, we set v\V2 = 1 and choose the 
parameterization 

Vl = ( U2 )-i = sin ( A ~^ ; p = 2 cos A. (2.20) 
sinw 

Then, using equations (2.13), (2.14) and 



i im ~^L=o = °> ( 2 - 21 ) 



/3-5-0 



one finds 



ahb=l 



1 ,. I W+-2/3 



lim [Z L 

/3->0 \ 



K(u) 



a^^a + ^ £ ^L a , 6 =o ) • (2- 22) 



This final result is quite different from equation (2.16) pertaining to ^ 0. This is not 
surprising since the criticality condition V\V2 = 1 imposed for ft = (Q = 0) is incompatible 
with the limit P — > (Q — )■ 0) of the similar condition uiu 2 = Q imposed for P ^ 0. 
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2.4 Enlarged periodic TL algebra 



Here we discuss the family of Enlarged Periodic Temperley-Lieb algebras £PTLn{cv, f3) la- 
belled by N and depending on the two parameters a = a^o and 0. Our discussion will show 
how the partition function in the loop model Zl can be calculated from the eigenvalues of 
the transfer matrix T{u) in certain link-state representations of SPTL^(a, 



Connectivities and the EPTLn(ol, (3) algebra: On a rectangle, N equidistant nodes 
are drawn on both the upper and lower edges. The rectangle is viewed as the planar repre- 
sentation of a band around a vertical cylinder, so that the left and right edges of the rectangle 
are identified along the so-called virtual boundary. A connectivity is then a connection of 
the 2N nodes by non-intersecting loop segments living between the upper and lower edges, 
such that every node is connected to one and only one other node. Loop segments can go 
through the virtual boundary, i.e. around the back of the cylinder. For instance, 

(2.23) 



is a connectivity. The product of connectivities is defined as follows: C2C1 is the connectivity 
obtained by drawing c\ above c 2 , reading the connectivity between the top and bottom 
N nodes of the ensuing connectivity diagram, and multiplying this diagram by the factor 
a Ua (3 ni3 , where np is the number of contractible loops and n a the number of non-contractible 
loops with homotopy {1,0} in the diagram formed by drawing c\ above c 2 . Because non- 
contractible loops do not appear if N is odd, the parameter a only comes into play for N 
even. As an illustration, we consider 



c 2 ci 




af5 



J \J \J \J V 



(2.24) 



The algebra SPTL^(a, (3) is the infinite-dimensional algebra formed by the linear span of 
the connectivities endowed with the product just defined. It is generated by the connectivities 
J, ei, e 2 , . . . , 6jv, Q and Q^ 1 : 



e N 
























n 



W—771 «"-jKZ3Si 



(2.25) 



where the index % on the TL generators is understood to be modulo N, with eo 



e N . 
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Products of connectivities then follow from the relations 



n 2 


= pe i} 












= ei_i, 




= ir 1 ^ = i, 








= e N , 


EQ ±X E 


= aE, 



> 1, 



(2.26) 



£ = e 2 e 4 . . . e N ^ 2 e N , 

where the last equation is for N even only. It is noted that £PTL N (a, (3) is actually generated 
by Q, Q^ 1 and any one of the TL generators e^. The enlarged periodic TL algebra is the 
quotient of the affine Temperley-Lieb algebra [10, 11,35,36] by the last relation in (2.26). 



The transfer matrix: The loop transfer matrix is an element of SPTL^(a, /3) given 
diagrammatically by 



N 



T{u)= - 


u 


u 






u 




> - 


> - 






^ - 



where 



u 



sin(A — u) 



+ sinw 



A— it 



(2.27) 
(2.28) 



— 2 cos A, and u is the spectral parameter. 



Link states and twist representations uja- Let iV equidistant nodes be drawn on a 
closed horizontal curve wrapping a cylinder. A link state is then a set of non-intersecting 
loop and line segments, drawn on the cylinder above the closed horizontal curve, linking the 
nodes pairwise or attaching vertical line segments to nodes. These vertical line segments are 
called defects and a loop segment between a pair of nodes is not allowed to connect above 
a defect. The number of defects of a link state is denoted by d, an integer that lies in the 
range 0, . . . , N and subject to the constraint N — d = mod 2. For instance, the link state 

w= -\ 1 1 r (2.29) 

has 2 defects. As indicated, it is convenient to depict link states on a horizontal line segment 
with the left and right ends identified. The set of link states with iV nodes and d defects is 
denoted by Bf^ and its linear span, Vfi, has dimension 

dimV^=^J. (2.30) 
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The counting of link states stems from a bijection previously exploited in [12] between link 
states with d defects and spin configurations (xi, . . . ,xjy) with G ±1 and J2i x i = If 
position i is occupied by a defect or if it marks the beginning of a half-arc (it connects to a 
position j by going towards the right), xi = +1. Otherwise i marks the end of a half-arc (it 
connects to a position j by going towards the left) and x% = — 1. The number of link states 
with half-arcs is equal to the number of states (xi, . . . , x^) with down arrows, and 
(2.30) follows readily. 

To define link-state representations of SPTL^la, /3), one must define an action cw of 
connectivities on link states that respects c 2 (ciw) = (c 2 ci)w for all connectivities Ci,c 2 and 
link states w. The action defining the twist representation u>d [12] is the following. First, the 
link state w is drawn above c, with the nodes of w connected to the TV nodes of the top edge 
of c. If any two defects of w are connected in the diagram, the result is zero. Otherwise, 
cw is equal to the link state obtained by reading the connection of the lower N nodes of c, 
multiplied by the scalar factor a na (3 ni3 v A where np and n a are respectively the number of 
loops in the diagram with homotopy {0} and {1, 0}, v is the winding parameter and A is the 
total winding of the defects, calculated as follows. If cw ^ 0, any defect of w travels across 
the connectivity and connects with a node on the lower edge of c. Let the defects of w be 
labelled by i, 1 < i < d. One can then calculate the distance A, G Z traveled by the defect 

A, = [initial position of i] — [final position of i\ , (2-31) 

where the nodes of the link state are in positions 1, . . . , N. In order for A, to measure a 
distance, we let the final position of i exit the interval 1, . . . , N if the virtual boundary is 
crossed. Defects traveling towards the left or right result respectively in positive and negative 
powers of v. The total winding is then the sum of the windings of all the defects, 

A = ^A,. (2.32) 

i 

As the number of defects is conserved, this defines a representation ui d of £PTL N (a, (3) for 
every allowed value of d. The non-contractible loop fugacity a comes into play in the Uq 
representation only, while the winding parameter v appears in all the representations with 
d > 0. Here are three examples of the action of cw. 

/3V 1 1 r\ (2.33) 



C\W 



.... Ij^S - 



(2.34) 



c 2 w 




a 



(2.35) 
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In the first example, the first defect has Ai = 4 and the second A 2 = 2, yielding a total wind- 
ing of A = 6. We stress that the winding parameter v is a parameter of the representation 
Ud and not a parameter of the lattice model. 

Closing the cylinder to form the geometric torus: Our interest in the enlarged 
periodic TL algebra SPTL^^a, ft) stems from the torus loop model discussed in Section 2.1. 
Indeed, we now proceed to show how the partition function of the loop model on the torus 
can be expressed in terms of T M \u), where T M \u) is defined as the (vertical) concatenation 
of M copies of the transfer tangle T{u). 

As an element of £PTLj^(a, ft), T M {u) is a weighted sum over all 2 MN configurations a 
appearing in Z^. Each such weight is given by an expression of the form 

a n "p n f>($mu) nm (sin(A - u)) nm (2.36) 

where n a and rip are the numbers of loops with homotopy {0} and {1, 0} that do not cross 
the upper and lower horizontal edges. 

However, this is different from the weight in of the corresponding configuration on the 
torus: Along with the factors of ft and a for the {0} and {1,0} loops that do touch the upper 
and lower horizontal edges, the factors a a f> for loops with homotopy {a, b} are also missing. 
(Here and in the following, we must exclude (a, b) = (1,0) in order to avoid counting twice 
the loops of homotopy {1, 0}.) To assign the correct weight to a configuration on the torus, 
we thus define the linear functional 

J= : £PTL N (a,ft) -> C, JF{c) = a n «ft n e J{ o>' b (2.37) 

a A 6=1 
(0,6)^(1,0) 

where an n counts the loops crossing the horizontal edge along which the cylinder is glued 
into a torus. It is evident that every loop with homotopy {a, b} ^ {1,0} crosses this edge, 
so n a>b = n^b- 
As an illustration of the evaluation of J 7 , let us first view the particular configuration in 
Figure 1 as a product of connectivities in SPTL^la, ft). The loop numbers in (2.36) are 
then n a = and rip = 2. Removing the rip = 2 contractible loops from the configuration 
yields the connectivity 




(2.38) 



and it readily follows that J-(c) = ft 2 a\ x . 

In (2.37), the linear functional T is defined algebraically as an operator mapping con- 
nectivities c G £PTLj^(a, ft) onto C. By construction, J 7 can be interpreted geometrically 
as gluing together the horizontal edges of a cylinder on which the loop configuration corre- 
sponding to a connectivity c is realized. Indeed, T has been designed such that the value 
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J~{c) € C is identical to the weight of the loop configuration on the torus resulting from this 
gluing process. 

It is recalled that the product cc' of two connectivities is itself a connectivity multiplied 
by a constant. Let us denote the resulting connectivity by c" and the constant by k. We 
thus have cc' = kc" and Tied) = kJ-(c") is the weight of the loop configuration on the torus 
corresponding to the product cc'. The constant k is given by the product of the fugacities 
of the loops in this configuration that do not intersect the horizontal edges of the cylinder, 
while the fugacities of the loops that do intersect the edges are incorporated in T(c"). Since 
e'e may differ from cc' = kc", the separation into the similar contributions coming from e'e 
may be different. On the other hand, because the two products cc' and e'e are equivalent 
as loop configurations on the torus, it readily follows that the two products cc' and e'e are 
mapped to the same product of fugacities by J 7 . We have thus established that J 7 has the 
cyclicity property 

J-(cc') = J-(c'c), c, c G SPTL N (a, f3). (2.39) 

Being also linear, T is therefore a trace operator on SPTL^(a, 0). 

By construction, the partition function of the loop model on the torus is given by 

Z L = F(T M (u)). (2.40) 

As we will see below, Z L can be expressed as a modified matrix trace in the twist represen- 
tations io d {T M {u)). 

2.5 Traces in the twist representations 

Our next goal is to write the linear functional J 7 in terms of traces in the twist representations. 
The main idea is to use the winding parameter v to decompose the traces of ujd{c) into Fourier 
modes and to assign the correct weights to each mode to produce J-"(c). Because T is linear, 
it can be applied to linear combinations of connectivities in £PTLn(c( } (3) such as T M (u). As 
written in (2.42) below, the expression for J-"(c) applies to connectivities c with a maximum 
of M strands crossing the virtual boundary. Here we write the Fourier decomposition in 
integral form and introduce the convenient notation 

Tr d (c) = Tru d (c). (2.41) 

The final result is the following proposition which we prove in the remainder of this section. 
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Proposition 2.1 The partition function for the loop model on an M x N torus is given by 



F( c ) = 7r dfj,y2g d (fj,,a)Tr d (c), 
27r Jo , 



Z L = F{T M (u)) where < 



(2.42) 



Go((i>,ac) =t = l, 

M 

Q d {n,a)= v~ Nk C k>d , (d>0) 

k=-M 

Ck,d = tkf\d{oi_k d), 

kAd 1 kAd 

I t n (a) = (2-5„, )T n (f). 

27ie winding parameter is parameterized as v = e lM anc? T n (x) is the n-th Chebyshev polyno- 
mial of the first kind, given by T n (cos8) = cos(n6). The role of the integral is only to select 
out the constant term in the Laurent expansion in v = e* 



3 *m 



For a given connectivity c, the matrix elements of ^(c) are powers of a and 0, but most 
of them are zero. In fact, in every column of Ud(c), at most one element is nonzero. Some 
matrices, such as wjv(ei), are simply equal to the zero matrix. To calculate Tr^c), we search 
for the link states in Bfj that are eigenvectors of c with eigenvalue different from (for 
generic a and (3). We will refer to such eigenstates as link eigenstates. 



An example: We start by working out an instructive example in detail: the connectivity 
c = n 4 e £PTL 8 (a,/3)i 

N=8 




As indicated in the illustration of this connectivity, four strands are crossing the virtual 
boundary. Thus, if Q 4 is realized on an M x iV torus (with N = 8), M necessarily satisfies 
M > 4, while the ensuing loop configuration on the torus has four loops with homotopy 
{1,2}. It follows that 

F(tf)\ N=s = aj )2 . (2.43) 

We now proceed to verify that the prescription (2.42) yields the same result. 

The link pattern with 8 defects is an eigenstate of c, with eigenvalue v 32 . A simple way to 
calculate the winding A of such an eigenstate is to count the number of loop segments con- 
nected to defects that cross the virtual boundary and multiply this number by N. Examples 
of other link eigenstates of c are ^11^11^" and f7\\ with eigenvalues v 16 and 1, 

respectively. Only a fraction of the four loop segments that cross the virtual boundary are 
connected to defects; ~ and for the first and second state, respectively. The connectivity 
c has no link eigenstates with 2 or 6 defects. For which sectors do these eigenstates exist 
and how many are there? If one takes any link state in Bf for d! = 0,2,4 and places two 
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copies of it side by side, thereby forming a state in B^ d ' , the result is an eigenstate of c with 
eigenvalue v 8d . Since this procedure to construct link eigenstates is exhaustive, the number 
of link eigenstates of Q 4 in the sectors d — 0,4 and 8 are simply the number of states in B®, 
B\ and Bf, respectively. It therefore follows that 

Tr 8 (fi 4 ) = v 32 , Tr 4 (fi 4 ) = Q v 16 , Tr (fi 4 ) = Q > Tr 6 (Q 4 ) = Tr 2 (fi 4 ) = (2.44) 

and subsequently 

JC ( q4 )\n=8 = ^^ 2 ^/i(^o(/^,a)Tr (fi 4 ) + ^ 4 (/i,a)Tr 4 (fi 4 ) + ^ 8 (/i,a)Tr 8 (fi 4 )). (2.45) 

The first term is easily evaluated to ( 4 ). In the second term, only the contribution for 
k = 2 in the sum in Q^fJ,, a) survives the integration, thereby yielding (f) 6*2,4 — 4i 2 (a:i )2 ) — 
4(af 2 — 2). In a similar fashion, we find the last term to be C 4j8 = £4(^1, 2) = a\ 2 — 4af 2 + 2 
and, overall, 

^(^ 4 )U =8 = 6 + 4« 2 - 2) + « 2 - 4a? i2 + 2) = < 2 (2.46) 

as required. 



For c a power of fl: If c = Q l for some % G Z (where Q° = J), the diagram of c has |i| 
loop segments that cross the virtual boundary. If tt 1 is a connectivity that contributes to Zl 
and T M (u), then M > When the top and bottom of c are connected into a torus, the 
number of loops is n = N A i and their homotopy is {a, b} = {i/n, N/n}, so that 

= a^ N/n , n = N A i. (2.47) 

The link state with only defects is an eigenstate with eigenvalue v m . In general, a link 
eigenstate of Q l is constructed by adjoining ^ copies of a link state with n nodes and r half- 
arcs (0 < r < [f J)- This gives |-B" _2r | = (") link eigenstates in the sector with d = ^(n — 2r) 
defects. Among the connections to the virtual boundary, only I! —^ L x \i\ are connected to 
defects. The eigenvalue is v Nl ( n - 2r )/ n an d 



n 




M("-2r)/n d=%(n-2r), r = 0,...,LfJ, 



Tr d (n t ) = {\rJ ' « v "'•••'L2J' (2.48) 

otherwise, 



while 



v Ni(n-2r)/n 



LtJ / x LfJ , . (2.49) 

= ( r ) %™-2r),£(n-2r) = ^ ( r J *n-2r(Oi/n,JV/n)- 
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For our prescription to work, this must reproduce (2.47), which is easily verified as a rela- 
tion between Chebyshev polynomials. In terms of trigonometric functions where oti/ n ^/n 
2 cos 0. the corresponding equality reads 



f n— 1 



(2cosfl) 



n 



= < 




n odd 



(2.50) 



n 



n even. 



This can be verified by an explicit computation of the Fourier series of (2cos#)™. 

The general case: We now verify that the expression for J-(c) given in equation (2.42) 
is correct for any connectivity c different from f2\ For such a connectivity, the state with 
only defects is not an eigenstate. Every connectivity has at least one link eigenstate, and 
among all the link eigenstates of c, one is particular: It has more defects than any other link 
eigenstate of c, and we will denote by j(c) the number of defects of this particular state. We 
now present an algorithm to determine j(c) for any connectivity c. Take for example the 
connectivity 



Because the nodes 4 and 5 on the bottom edge are connected by a loop segment, any link 
state w satisfying cw = Xw for some A G C will have the nodes 4 and 5 connected. The same 
argument applies to the nodes 7 and 8, and w = 111^1^, where we are using question 
marks "?" to indicate the nodes whose connections are not yet determined. When this w is 
drawn atop c, a factor (3 appears because a loop is closed at positions 7 and 8, and the nodes 
2 and 3 of the bottom edge are connected, so w = 1 n r\ 1 n . Now, cw is seen to produce 
a factor of /3 2 , and no further connections are made between nodes on the lower edge of c. 
The unfixed nodes that remain can be replaced by defects, and their number is j(c) = 2. 

The general procedure is as follows. If nodes on the lower edge of c are connected, the 
corresponding nodes must also be connected in any eigenstate w of c. Acting on the partially 
known w may then close new connections between nodes on the lower edge, which are in 
turn added to w. This process is repeated until no new connection is formed between nodes 
on the lower edge. The nodes of w which are left undetermined can be replaced by defects 
and their number is j(c). The resulting link state is clearly unique. 

A small miracle operates here. First, the link state obtained from this algorithm has 
eigenvalue a na (3 ni3 v A for some n a , ftp and A, where A depends on the permutation of the 
unfixed nodes. Second, the factor a Ua fi nfi is precisely the weight of the loops with homotopy 
{0} and {1,0} that intersect the upper/lower edge of c when the torus is formed by gluing 
together the horizontal edges of the cylinder. 

If the above procedure for the connectivity c fixes every node of w, then Tr^(c) = 
for d > 0, w is the unique link eigenstate of c and Tr (c) = a na (3 ni3 . It subsequently 




(2.51) 
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follows that c has no loops winding around the torus vertically and from equation (2.42), 
JT( C ) = J o 2 ^Tr (c) = a na (3 n ? as expected. 

If a subset of the nodes of w are undetermined, the positions of these nodes are those 
in c where loops with homotopy {a, b} cross the horizontal edges. Up to a deformation of 
the loop segments, the algorithm produces a diagram where the connectivity c acts on the 
j(c) nodes with question marks by operating a cyclic permutation, as Q l did in the previous 
example. For example, for a given connectivity (different from (2.51)), this diagram could 
be 




Replacing the question marks of w by defects gives the unique link eigenstate with j(c) 
defects. The shift of the j(c) nodes is denoted by i(c) and is measured by the number of loop 
segments that cross the left /right virtual boundary, \i(c)\. On the torus, the connections 
are divided into n = i(c) A j(c) loops of homotopy {i(c)/n, j(c)/n}. If n > 1, c will have 
other link eigenstates, where some of the question marks are replaced by half-arcs instead of 
defects. The picture is quite similar to what we encountered for VL l : Separate the j(c) nodes 
of w into j(c)/n sequences of n nodes, choose a link state with r half-arcs (r = 0, 1, . . . , |_§ J) 
and adjoin j(c)/n copies of it side by side to form a link state with d = — (n — 2r). The 
result is an eigenstate of c, and because only the fraction lt ^ L of the \i(c)\ loop segments 
crossing the boundary are connected to defects, the associated eigenvalue is /3 n Pv Nl ( c ^ n ~ 2r ^ n . 



The number of such link eigenstates in the sector with d 



3(c) , 



n 



2r) defects is (™) and 



^ n2ir L 2 J 

— / d/j,y^g J{c) 

7T In L ' 



2 

LfJ 



(n-2r) 



(fj,,a) 



n \ v Ni(c)(n-2r)/n 



- 1 (: ! ) c H 

f3 np (0!i(c)/n,j(c)/n) 



^(n-2r),^i(n-2r) 



L"J 

p x ± (;) *„- 



2r(0!i(c)/n,j(c)/n) (2.52) 



where the last equality follows from (2.50). Because the expression for T(c) is linear in c, 
this completes the proof of Proposition 2.1. 



3 Critical dense polymers 

An objective of this work is to express the partition function of the loop model on the torus 
in terms of the transfer matrix, and to test whether the result is modular invariant. From 
now on, we focus on (3 = 0, the value corresponding to critical dense polymers, while we 
initially keep the parities of M and iV and the values of the non-contractible loop fugacities 
ccjj free. In this general setup, the transfer matrix satisfies an inversion relation [7] which we 
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solve for the twist representations w^. As shown in Section 2.5, the problem of computing 
the loop partition function boils down to calculating the eigenvalues of the transfer matrix in 
these representations. As a result, we obtain the general expression for the partition function 
of the loop model for arbitrary M, N and a^-. Considerable simplifications are subsequently 
possible if we set all the non-contractible loop fugacities to aty = 2. In this case, we prove 
that the partition function is modular invariant if both M and N are even. 

The random cluster model, which is defined here only for M and N both even, is also 
expected to yield a modular invariant for some choice of the free parameters. It is recalled 
that, although both the loop model (Section 2.1) and the random cluster model (Section 2.2) 
are well-defined for (3 = \J~Q = (and give nonzero partition functions), the passage from 
one to the other in Section 2.2 is awkward. Consequently, the FK partition function we 
calculate is of the no-cycle ( "restricted" ) random cluster model with Q = w + = discussed 
in Section 2.3. 

Of course, setting the weight of nontrivial clusters to zero gives a trivial modular invariant, 
Z' FK = 0. This is obviously of very little interest, so one of our goals is to build another, 
nontrivial one from the transfer matrix approach to the loop model. With (3 = 0, this 
method will allow us to compute the first term in (2.22), but not the second one since this 
would require knowledge of the spectrum of the transfer matrix away from /3 = 0. Evidently, 
a direct way to eliminate the second term in (2.22) is to set the weight of cross-topology 
clusters to zero, w + = 0, in which case the restricted FK partition function follows from that 
of the loop model. 

Let us comment on some of the terminology we will be using. The number of defects 
d classifying the representation ujd is a quantum number that separates the theory into 
Ramond (~ even), Neveu-Schwarz (| odd) and Z 4 (d odd) sectors. Following [7,19,37], we are 
thus using terminology of supersymmetry even though we do not claim any superconformal 
symmetry in our model. 



3.1 Spectrum of the transfer matrix 

For (3 = 0, it was shown in [7] that the transfer matrix satisfies the inversion identity 

T(u)T(u+§) = (cos 2N u + (-1) N sin 2N u) I + (cosusmu) N J 
where J is the vertical identity diagram while J £ £PTLj^(a, = 0) is given by 

N 





with 







1 

















(3.1) 



(3.2) 



Using the Drop-Down Lemma in [7], we find 

(_l)V((-l)V 7V + ^- 2JV ), d>0, 
2 2), 



U d (J) = UJd(I) x 



d = 0. 



(3.3) 
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From the inversion relation, exact expressions for the eigenvalues of the transfer matrix can 
be found. Note that if a is replaced by v N + v ~ N , the expression for 0Od(J) for d > also 
applies for d = 0. Therefore, all the results for d > 0, where the eigenvalues of 0Ud(T(u)), 
denoted by T(u,/i), depend on the winding parameter v = e* M , can be extended to d — if 
we set 

a = 2 cos fiN. (3.4) 
Then, in the sector with d defects, we have 



u d (T(u)T(u+ f)) = oj d (I) x F(u,v) 



(3.5) 



where 



F(u, v) = (v N cos N u + (-1)tV n sin N u) (v~ N cos N u + (-l)^v N sin N u) . (3.6) 

The expression for F(u, v) can be factorized, 

F(u,e i n=T(u,fi)T(u + i,ri, (3.7) 

and the possible values for the eigenvalues are thus found by separating the zeros into 
contributions from T(u, fx) and T(u + |, /i). For iV odd, we find 

F(u, e >n = e ' 2Msi ° 2( |;. ( r 1) ^^ ) n(^ + .«w + (-1^*4 <3. 8) 



i=i 



d+l 



T(u,fi) = 



ee-T e -» w sin(f + (-1) 2 ^iV) 



N 

n(e 2 - + ^tan(^ + (-l)^)) : 
j'=i 



(3.9) 



where e, ej G {+1, —1} for all j — 1, . . . , TV. For TV even, we find 



< e~ 2iuN sin 2 (i^N) S 



F(u, e*) = < 



2 2N-2 n( e4m + tan2 (f -»))> 

^n(e- + tan 2 ( 



2 { (?]-!)■* _ .. 
2N V" 



T(u,fi)={ 



e(-i)2e~ iuN sin(//TV) 
2^i 

e(-i)^e~ iuN cos(/zTV) 

OJV-I 



J'=l 



no 



e 2m + iej tan 



j'=i 
v 



f odd, 
| even, 



| odd, 



(3.10) 



(3.11) 



n(e 2 - + ^tan(^-/,)) 
i=i 



\ even, 



again with e, 6j G {+1, —1} for all j — 1, . . . , TV. For both parities of TV, the eigenvalues are 
characterized by the choices of e and €j, j = 1, . . . , TV. This gives a total of 2 N+1 possible 
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eigenvalues. Because a given sector only has y^y eigenvalues, only a subset of the possible 

eigenvalues appear in the sector d. This information is encoded in the selection rules of 
Section 3.2. 

The eigenvalues found in [7] are for \i = 0, and extra care has to be taken in the Neveu- 
Schwarz sector (N even and | odd). Naively, the result for T(u,/j) in (3.11) appears to be 
zero in that case because of the prefactor sin (/iiV), but the j ' = y term of the product is 
singular in the limit /i — > 0. The correct expression is obtained by applying l'Hopital's rule. 
More generally, similar resolutions apply if \x is such that sin(^ + (— l)~[j,N), sm(fiN) or 
cos(/xiV) is zero. 

To describe all three sets of eigenvalues in (3.9) and (3.11), we will henceforth use the 
short-hand notation 

N 

T(u,n) = K(fj) x e x Y[ (e iu + ie- iu e j tanx j ^, (3.12) 

i=i 

where 



inN 

e ' * 



sin(f + (-l)*J>JV), iVodd, 



K(ji) = x I sin/iAT, N even, § odd, (3.13) 
cos /jlN, N even, | even, 

fi^ll + AT odd, 

Sj- = {i£ - /x, AT even, \ odd, (3.14) 
(23 2 ^ )7r - /i, A^ even, | even. 
It is noted that the eigenvalues (3.12) satisfy the crossing symmetry 



T(Z-u,fi)=T(u,ti), (3.15) 
here indicated for complex u and \x. 



3.2 Selection rules 

In this section, we write down the selection rules dictating which of the solutions to the 
inversion identity appear as eigenvalues in the spectrum of Ud(T(u)). The multiplicities of 
the eigenvalues also follow from the selection rules. In [7], the corresponding selection rules 
for /i = are expressed in terms of column configurations and excess parameters a and a. 
Here they are formulated exclusively in terms of conditions on the parameters ej and e, and 
their proof is given in Appendix A. That the selection rules reduce correctly to those of [7] 
for jj, = is also verified in the appendix. 

Another component of the selection rules concerns the identification of the groundstates. 
Eigenvalues of Ud(T(u)) are in general complex, as can be seen from T{u — 0) = Q whose 
eigenvalues in the sector with d defects are of the form T(0,/i) = v d e^ ZL , for some m G Z. 
Hereafter, the groundstate eigenvalue of uJd{T(u)) is defined as the one with maximal norm 
and is denoted by Tj. In Appendix A, Td will be identified as the eigenvalue for which the 
linear term, Eig("H), is maximal. 
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Z 4 sector (N odd): 

(i) In the sector with d defects, the parameters ej, e G { — 1, +1}, j = 1, ■ ■ ■ , N, character- 
izing an eigenvalue of the transfer matrix satisfy 

N 

X)(-l) i e i = d(-l)^ > 6=(-l)^. (3.16) 

The number of different choices for the set of parameters {e 3 -} satisfying (3.16) is 
precisely (il^) , and each appears exactly once. 

(ii) Let 6j = tN+i-j for j = 1, . . . , ^j^. For \i G j^}, the groundstate in the sector 
with d defects is characterized by €n+i = +1 and 

I ^ A > J 2 ' • • • ' 2 ' / 01 m 

If we set z/fc = (— l) fe+ ~^~efc, the selection rules simply state that exactly of these 
N renormalized parameters satisfy V\. = —1. The full spectrum of Ud(T AI (u)) can thus be 
written using a generating function as 

m , M , xN 1 fd n G d (z,v)\ N-d , 

MT M (u)) = - I AJ L ) n = __ (3.18) 



2=0 



where u = and 

N 



G d (z,v) = {K(ji) x e) M Y[ ((e iu + ie~ iu tanx' j ) M + z(e iu -le-^tanx'j)^ , (3.19) 



3=1 

(3.20) 

Ramond and Neveu-Schwarz sectors (iV even): 

(i) In the sector with d defects, the parameters ej, e G { — 1, +1}, j = 1, • • • , N, character- 
izing an eigenvalue of the transfer matrix satisfy 

N 

5> = -d, C =(-1)W. (3.21) 

3=1 

The number of different choices for the set of parameters {e.,} satisfying (3.21) is 
precisely (e-aJ , and each appears exactly once. 
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ii) For \i G the groundstate is characterized by 



1 i — I i I -+- 1 iV I d+2 I 

, - a, J LiJ + V--, 2 4 (3.22) 
— 1, otherwise. 



The full spectrum of oo d (T M (u)) admits the expression 



M T- { n)) = I n= ^A, (3.23) 

1 ' n\ V dz n J z=o 2 ' v ; 



where v = e* M and the generating functions are defined by 

N 



G d (z,v) = (K(fi) x e) M \\ ({e iu - ie~ iu taxi Xj ) M + z (e iu + ie~ iu tanx,) A/ ) (3.24) 



with Xj given in (3.14). Note that this expression differs from the similar one in (3.19). For 
later convenience, it is also noted that 

G 2 (-l,l) = limG 2 (-l,e*) = 0. (3.25) 
3.3 Lattice partition function 

We now proceed to find an exact expression for the partition function of the loop model 
from the generating functions G d (z,v). First, we make two crucial observations: 

(i) As a straightforward calculation shows, 

G d (z,e^ + ^) = (-l) M G d (z,e^). (3.26) 

Indeed, K(fi+ |U = —K(fi) and the terms in the product (3.19) (for N odd) or (3.24) 
(for N even) are mapped onto one another when [i is shifted by jr. 

(ii) G d (z, v) admits a Laurent series expansion in the winding parameter v with finitely 
many nonzero coefficients, with the terms with minimal and maximal powers of v given 
by %r NM and v NM , respectively. 

Indeed, for N odd, for example, the generating function in (3.19) can be expressed as 

N 

G d {z,v) = (e"^e) M ]^ (^(e^ cosM + 2sinue-^) M + z(e-^ cosM + zsinue^) M ) , 

3=1 

(3.27) 

where (A. 15) has been used to simplify the factor K M {fi). Because 

& = ^xp(^(2j " 1)(-1) J+ ^), (3-28) 
the result readily follows. Using (A. 27), a similar argument carries through for iV even. 
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It follows from these two observations that G d (z,v) is v NM times a Taylor series in v 2N . 
Because G d (z,v) is also polynomial in z, it can be written as 



M N 



G d (z,v)= Yl X>*Vft,(*,*). (3.29) 

k=-M,-M+2,... 1=0 

In terms of the coefficients g d (k,£), the trace of u d (T M (u)) is 

M 

Tr d (T M (u))= v Nk g d (k,^). (3.30) 



k=-M-M+2,. 



One can now use this last equation along with equation (2.42) to write the partition function 
of the loop model as 



N M 



E E Ck,d9d(k,Z=±), N odd, 

d=l,3,... k=-M,-M+2,... 

N M 

Tr (T M ( M ))+ E C kid9d (k,^), N even. 



d=2,4,... k=-M-M+2,... 

It is observed that the integers k and M always have the same parity. This is not a surprise 
since, in Section 2.5, k counts the number of times defects cross the virtual boundary, and 
this number always has the same parity as the number M of layers of T(u). 

The special term Tr (T M (w)), which appears for N even, can be rewritten as follows. 
First, it is recalled from (3.4) in Section 3.1 that the parameter a = «i,o is conveniently 
parameterized as v N + v~ N . Because Tiq(T m (u)) is polynomial in a, the coefficients in 
equation (3.30) for d = must satisfy the symmetry 

g (h,%)=g (-k,%). (3.32) 

(This also follows readily from equation (3.52) below.) We thus have 

Tr (T M (n)) = | £ (2 - 5 k>0 ) (v Nk + v~ Nk ) g (k, f ) = J] C k , Q g (k,f), (3.33) 

0<fc<M 0<fc<M 
fc=Mmod2 fc=Mmod2 

where C fe , , < k < M, is defined as C M for d > in (2.42): 

„ jt k (a lfi ) =t k (a), k>0, 

Ck ' = \i, k = o. (3 - 34) 

Equations (3.31) and (3.33), combined with (3.19) or (3.24) for N odd or even respectively, 
give the general expression for the partition function of the loop model. 
Considerable simplifications are possible if we set 

atij = a = 2 (3.35) 
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for all In this case, t n = 2 — 8 n>0 and, from equations (3.31) and (3.30) at v = 1, 

^L,= 2 = E (2-S d , )Tr d (T M (u))\ a=2 . (3.36) 

0<d<N v=l 
d=N mod 2 

For each parity of the system size N, the total number of contributing link states is thus 
given by 

(2-5 dj0 )dimV$ = 2 N (3.37) 

0<d<N 
d=N mod 2 

where the dimensions are given by the binomial coefficients (^-^j as in (2.30). The explicit 

expression (3.36) giving the partition function as a sum of ordinary sector traces is key to 
the identification of the corresponding modular invariant discussed in Section 3.5. 

3.4 Finitized partition functions 

The results of Section 3.3 are expressed in terms of the N + 1 functions G d (z, v) where 
d = 0, . . . , N. In the following, we show how to extract the physical information from the 
three functions G (z,v), Gi(z,v) and G 2 (z, v) alone, and find that this allows the direct 
evaluation of the finitized partition function for M even and a^j = 2 for all i,j. In contrast, 
similar results have not been obtained for M odd. It is noted that non-contractible loop 
fugacities are also set equal to 2 in the work [19]. 

The expansion of the generating function G d (z,v) as a Taylor series in the variable z is 
given by 

N M 

G d (z,v)=J2 A *} zi > A S= E vNk 9 d (k,£) (3.38) 

£=0 k=-M,-M+2,... 

and is such that 

Tr d (T M (u))=A^ (3.39) 

a, 2 

Recalling that T d denotes the groundstate eigenvalue of u d (T(u)), we introduce normalized 
generating functions, 

G d (z,v) = ^^\, v = e» (3.40) 

Now, to compute the finitized partition function, we first determine the appropriate 
finitizations G d N \u, v) of the normalized generating functions. This requires restricting fi to 
be of the order 1/N, and we choose to parameterize it as 

A*=^, aeH.il, (3-41) 

as in the description of the selection rules for the groundstates in Section 3.2. The finitizations 
G d N \u,v) are then obtained following the approach of [7] based on the evaluation of finite 
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excitations. The results for d = 0,1,2 are discussed below and are expressed in terms of the 
aspect ratio 

M 

I = ¥ (3-42) 
and the modular nome q and its complex conjugate q, given by 

q = exp(-2vr^e- 2iu ), q = exp{2m5e 2iu ). (3.43) 

Here it is recalled that the models we study are defined on a rectangular array of Mx JV 
square tiles with periodic boundary conditions to cover the torus. In the continuum scaling 
limit, the periodicity of the torus is described by the modular parameter r which lives in the 
upper half of the complex plane, i.e. Im(r) > 0. The partition functions we compute are 
expressed in terms of a modular parameter r through the associated nome q and its complex 
conjugate q defined by 

q = e 27TiT , q = e- 27Tif . (3.44) 
Comparing this with (3.43), we see that 

r = -Se' 2tu . (3.45) 

As discussed at the end of Appendix B, the spectral parameter u (for < u < ~) can be inter- 
preted [38] geometrically as a measure of the spatial anisotropy of the lattice thus explaining 
its appearance in the definition of the modular parameter. For isotropic interactions, u — | 
and 

T = -5e-™' 2 = % ^ (3.46) 

as expected. 



Z 4 sector: For N odd, the generating functions have the symmetries 

(n e N), G 3 (z, v) = z N (-l) M G 1 (z- 1 , v' 1 ). (3.47) 



G 4n+1 (z,v) = (-l) nM G 1 (z,v), 
G 4n+3 (z,v) = (-l) nM G 3 (z,v), 



This allows us to write 



I^M^W d-l 
1J 4 A N - d , ~2~ 



even 



Tr d (T M (u)) = { ^ ,\\ (3.48) 

The partition function for the loop model with M even and ctij = 2 (see equation (3.36)) 
then simplifies to 

N 

d=l,5,9,... " d=3,7,ll,... '" e=0 

= G 1 (l,l) + (-l)^G 1 (-l,l). (3.49) 
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As outlined above, with // as in (3.41), we work out the finitized generating function 

3=1,3,... j 2.1.... 

x J] (l + z{-l) M q^- a ) J] ((-l) M q^ +a + z) (3.50) 

3=1,3,... j=2,4,... 

where the modular nome is defined as in (3.43). For M even and iV odd, we then obtain the 
finitized partition function 

N+l N-l N+l N-l 

Z { L N) \ aiJ =i = ^ox(f[(l + ^) n(l + ^)+fl(l-^) fl(l- f^))- (3.51) 

3=1 3=1 3 = 1 3=1 

The number of contributing link states is obtained by setting q — q — 1 and is readily found 
to be 2 N as in (3.37). 

Ramond and Neveu-Schwarz sectors: For iV even, the generating functions have the 
symmetries 

G An (z, v) = (-ir M G (z, v), G (z, v) = z N G (z-\ tr 1 ), 

G in+2 (z,v) = (-ir M G 2 (z lV ), [H h G 2 (z,v) = z N (-l) M G 2 (z-\v- 1 ) K ' } 

from which it follows that 

-l)i M A^, I even, 

Tr d (T M (u)) = { \\ ' (3.53) 

\-l)^ M A^, | odd. 

> 2 

Further simplifications are possible for v — 1 and M even as the second set of relations in 
(3.52) then translates into A^j = A^ N _ l for i — 0, 2. It subsequently follows that 

JV N_ 
2 2 

G t (l, 1) = £(2 - )4 1 £ ) , Gtf-l, 1) = £(2 - 6 e » ){-lfAf} (3.54) 

£=0 £=0 

for i = 0, 2, and the partition function for the loop model can be expressed as 



d=0,4,8,... " d=2,6,10,... 

= i (g (1, 1) + (-l)TGo(-l, 1) + G 2 (l, 1) - (-I)t G 2 (-l, 1) ) (3.55) 
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where the under-braced term vanishes due to (3.25). With the parameterization (3.41), 
the fmitizations of the (normalized versions (3.40) of the) generating functions appearing in 
(3.55) are given by 

Gf\z,e^) = fj {l^~ a + *) (1 + zq^ +a ) 
i=i 

LfJ 

x fl{(-l) M q M ^ +a + z)(l + (-l) M zq 2J ^- a ), (3.56) 
i=i 

LfJ 

G { 2 N \z,e^) = (l + zq a )(l + (-l) M zf a )f[(q j - a + z)(l + zq j+a ) 

3=1 

L^J 

x f[ ((-l)V'+ a + ^)(l + (-l) A W" a ), (3.57) 

3=1 

where the modular nome is defined as in (3.43). Using the result 

(?# (3.58) 

obtained in [7] from an Euler-Maclaurin analysis of the eigenvalues for p, = 0, we can now 
write the finitized partition function for M and iV even as 

4a, jrf =Trux (in (i+^rKi+f^) 2 

i=i i=i 

I N+2 j |^JVj |^JVj Y N ~ 2 \ 

+ ltl (l-^) 2 f[(l-^) 2 + %?) l ri(l + ? i ) 2 II (1 + ^) 2 )- (3.59) 

j=l 3=1 j=l 3=1 

The number of contributing link states is obtained by setting q — q — 1 and is readily found 
to be 2 N as in equation (3.37). 

An interesting parallel can be drawn between the calculations of the finitized partition 
function of critical dense polymers and that of the Ising model carried out in [39]. The spin 
transfer matrix of the Ising model also satisfies an inversion identity from which its spectrum 
can be determined. The resulting finitized partition function of the Ising model is of a form 
very similar to the one of critical dense polymers (3.59) and also contains three terms. In the 
context of Q-state Potts models, the inversion identity satisfied by the spin transfer matrix 
is a feature unique to the Ising model (Q = 2). Likewise in the context of the logarithmic 
minimal models CM.{p,p'), we believe that the loop transfer matrix satisfies a functional 
relation of degree p', implying that the inversion identity (3.1) is unique to critical dense 
polymers CM.(1, 2). We hope to discuss this elsewhere. 



lim 



rpM 

1 2 



M,N>1 T^ 1 
M=8N 



fi=0 
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The analysis of the finitized partition functions above is based on a torus formed by 
gluing together a rectangle without twisting the boundaries. We generalize this to helical 
tori in Appendix B where we find that, for ctij = 2 and even helicity t, the finitized (twisted) 
partition function is given by the expression (3.51) or (3.59), depending on the parities of 
M and N, but with q replaced by the nome in (B.9). 

3.5 Modular invariant partition function 

It is quite remarkable that the partition function (3.36) of the loop model on the torus for 
which all the non-contractible loop fugacities a^j are set to 2 can be written exclusively 
in terms of traces in representations with v — 1. In [7], the behaviour of these traces was 
studied in the M, N — > oo limit. Here we demonstrate that the partition function (3.36) is 
modular invariant in the continuum scaling limit if M and N are both even. It is recalled 
that the modular transformations are generated by 

S: r ^ , T : + 1 (3.60) 

r 

acting on the modular parameter r given in (3.45). 
Let us define Z d by 

TT d (T M (u))=T^{ x xe M xZ d , (3.61) 

where T max is the maximal eigenvalue of T(u), that is, T max = T for N even and T max = T\ 
for iV odd. In [7], it was found that these partition functions correspond to conformal 
partition functions associated with a logarithmic CFT with central charge 

c=-2 (3.62) 

and conformal weights given in terms of the defect number d by 

t 2 - 1 d / \ 

A* = -g-, t = -. (3.63) 

The corresponding Kac labels r, s, where t = 2r — s, are integer in the Ramond and Neveu- 
Schwarz sectors (N even) and half- integer in the Z 4 sector (N odd). 
We recall the specialized Jacobi theta functions 

*.M = E ^ ( 3 - 64 ) 

and the Dedekind eta function 

oo 

V(q) = q*T[(l-q n ), (3.65) 

n=l 

where the modular nome is defined in (3.43). For the simplified expressions in Section 3.4 to 
apply, we let M be even, in which case e M = 1 for all d. From the analysis of the continuum 
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scaling limit in [7], and up to the divergent term T^ x , we then find that the conformal part 
of the partition function (3.36) is given by 

_ itf l2 ( 9 )i 2 + i*» 2 ( 9 )i 2 

£ *(») = ' L p" . 

<fe2N-l 1 

for A" odd, while for iV even, it is given by 



Z (q) + 2 ZM = T-r-rr 2 £ \^M\ 2 - (3-67) 

d£2N ' 'W s=0 

As already observed in [7], the partition function (3.66) for A" odd is not modular invari- 
ant. This is not a surprise since the sectors with A" odd are incompatible geometrically with 
the sector containing the vacuum. The sesquilinear form in (3.67), on the other hand, is a 
well-known modular invariant Coulombic partition function [18,19]. The conformal partition 
function of critical dense polymers on the torus with M and A" both even is thus given by 
this modular invariant. Writing it as 

Z (q) + 2 Z d (q) = |x_i(<?)| 2 + 2|xo(g) + T + |Xf(g)| 2 , (3.68) 

dS2N 

where the W-irreducible characters [23] (see also [40]) are given by 

x_i(g) = ^ ^(g), xo(g) = ^{^M + v 3 (q)), 

i i 

%(g) = ^y^, 2 (g), xi(g) = ^(^1,2(9)-^)), 



(3.69) 



it is recognized as the modular invariant partition function of symplectic fermions on a Z2 
orbifold [21]. From Section 2.3, this is also the partition function of the no-cycle FK cluster 
model at Q = w + = in the continuum scaling limit. 

The W-characters incorporate the W-extended symmetry of the triplet model associated 
with symplectic fermions. As discussed in [41], the decomposition of the W-irreducible 
characters in terms of irreducible Virasoro characters reflects that an arbitrary number of 
defects, compatible with the system size, is allowed in the corresponding loop model. Indeed, 
this is what we find in the expression (3.68) for the modular invariant. 

We stress that the modular invariant partition function in (3.67) is not obtained by the 
usual matrix trace summing over all even A^ sectors since the gluing procedure governed by 
the linear functional T has changed the multiplicity of all of the d ^ sectors by a factor of 
2. Indeed, the standard matrix trace does not yield a modular invariant. In stark contrast, 
the modified matrix trace in (3.67) does yield a modular invariant partition function. 

4 Conclusion 

Because of the nonlocal nature of the degrees of freedom of the logarithmic minimal models 
£M(p,p'), the structure of their bulk logarithmic CFTs is much more complicated than in 
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the case of rational theories. In particular, the bulk representations are no longer given as 
simple tensor products over the left- and right-chiral halves of the theory. Similar observa- 
tions apply [42] in the context of the fj[(l|l) super-spin chains. In the case of critical dense 
polymers CM.(1,2), the bulk CFT of the associated Wi,2 (or ci^) triplet model was con- 
structed from the boundary CFT in [28] and, for the case of critical percolation £Ai(2,3), 
the bulk CFT of the associated triplet model was constructed from the boundary CFT 
in [43]. Generalizing the results of these papers, the modular invariant partition functions 
of the £Ai(p,p') models were conjectured in [44] from a lattice perspective. Although the 
representations do not factor, the modular invariant partition functions do factor in terms of 
characters and are given [44] as a combination of Coulombic and (rational) minimal modular 
invariant partition functions. To understand the bulk theory better, it is desirable to have 
a lattice derivation of these modular invariant partition functions. 

In this paper, we have taken the first step in this program by calculating, via an exact lat- 
tice derivation, the modular invariant partion function of critical dense polymers £Ai(l,2). 
Starting with the work of Cardy, there are well established physical principles which assert 
that, with purely periodic boundary conditions (and no seams or mismatched size parities), 
the torus partition functions should be modular invariant. These arguments apply equally 
to rational and logarithmic CFTs. Indeed, the modular invariance of the torus partition 
function of critical dense polymers is confirmed by our exact lattice calculation. In the lat- 
tice setting, the difficulty of closing the cylinder to the proper geometric torus is overcome 
by introducing a modified trace. In the continuum scaling limit, we obtain the expected 
modular invariant partition function of symplectic fermions on a Z 2 orbifold. The modi- 
fied trace acts to combine the three W-sectors from the cylinder with integer degeneracies 
(1,2, 1) as in (3.68). These degeneracies coincide [44] with Coxeter exponents given by the 

(2) 

right eigenvector of the Cartan matrix of the twisted affine Dynkin diagram A 2 with three 
nodes. The appearance of non-simply laced affine graphs and such degeneracies seems to be 
a feature of logarithmic CFTs. As is made clear by the lattice derivation, it is precisely the 
requirement of implementing a modified trace that leads to the correct degeneracies. From 
within CFT itself, and in the context of the general £Ai(p,p') models, it is still unclear 
precisely what principles select out the correct modular invariant partition functions. 

Our modified trace (2.42) holds for all of the logarithmic minimal models CM(p,p'). The 
problem of calculating the modular invariant partition function of the general £M(p,p') 
model is thus reduced to the problem of calculating the eigenvalues of the transfer matrix 
on the cylinder in sectors with d defects. We hope to come back to this problem in later 
papers. 

Helical boundary conditions can be applied to the construction of the torus by inserting 
a horizontal seam consisting of a power of the shift operator Q, as discussed in Appendix B. 
Moving beyond the modular invariant torus partition function, it is possible to introduce 
vertical integrable (topological) seams yielding general toroidal conformal partition functions. 
For rational theories, the fusion of these vertical seams or their corresponding twist operators 
is described [45] by the Ocneanu algebra [46]. It is an interesting question as to what extent 
the Ocneanu algebra needs to be generalized in the logarithmic setting, but this question is 
beyond the scope of this paper. 
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A Proof of the selection rules 

Here we prove the selection rules in Section 3.2. For N = 1, the transfer matrix is one- 
dimensional and its single eigenvalue is given by 

T(u, jj) = v cosm + v ~ 1 sinw. (A.l) 

This is readily seen to be reproduced by the selection rules which state that e = e\ = 1. To 
prove the selection rules for N > 1, we will study the first terms in the Taylor series of T(u) 
around u = and compare their spectra with the similar objects in the periodic XX model. 
The expansion of T(u) around u = is 

N 

T{u) = n(l + uU + 0{u 2 )), n = J2 e i> ( A - 2 ) 

1=1 

where we are using the sign convention for % of [12] for ease of comparison. By expanding 
T(u, n) around u — 0, one finds the eigenvalues of u^iT-L) to be 

( ^ + (-1)^, N odd, 
Eig(H) = e i sin2a ^ x i = \ f - A*, N even, f odd, (A.3) 

I - 2N ' fj,, JM even, | even. 

Indeed, 

N , N 



T(u, a) = K(u) x e x TT TT ( e cos Xj + ic e 7 - sin Xj) 

- LJ - cos Xh - LJ - V / 

k=l k j=l 

N ietxi N 

TT ( cos m + i sin ue' 2iejXj ) , (A.4) 

cosxl. - LJ - V / 



, . cosx k . 

k=l 3=1 



Eig(n) 
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where the under-braced prefactor is the corresponding eigenvalue of ^(fi). The eigenvalues 
of UdiH) are then 

N N 

Eig(H) = i ^^(cos 2xj — i€j sin 2xj) = ej sin 2xj, (A. 5) 

because the first term sums to zero with the Xj values given in (A. 3). As shown in [12], the 
spectrum of 0Jd(H) (for a = 2 cos fiN and — 0) is the same as the spectrum of the periodic 
XX Hamiltonian 

N 

H = J2 (^K+i + e- 2 ^++i) (A.6) 

in the sector with S z = |, where 

v 

5 *=lX>l> ^+1 = ^1. o €{+,-,«}. (A.7) 

To prove the selection rules, we proceed as in [5] by comparing the parameterizations of the 
spectra of the two hamiltonians H and H. First, we apply a Jordan- Wigner transformation 
to the XX Hamiltonian H. 

Proposition A.l In the sector with S z = |, H can be realized in terms of fermionic 
operators as 

F = ^2cos( 9 + 2^)77^ (A.8) 
q 

where 

( (2m+iv ^+d e ven 
1W = \ i ' ^ ' m = 0,...,N-l. (A.9) 

^ N ' 2 uuu ) 

The fermionic operators rj q and obey the usual anti- commutation relations 

0?J> Vq'} = <W, {Vq, Vq'} = {rfq, v\>] = 0, (A. 10) 

and are given in terms of Pauli matrices by 

N N 

3 q VN^ 3 

(All) 

k=l k=l 



The proof for ji = appears in [47]. Its generalization to the present case is straightforward 
and therefore omitted. 
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Using J2 q cos(g + 2/i) = 0, we immediately write 



# = -^2cos( g + 2/i)v7t. (A.12) 
q 

The state \U) = | + + + . . . +) is readily seen to be an eigenstate of H with eigenvalue 
0. More generally, a basis of eigenvectors in the sector S z = | is generated by the states 
Vqi ■ ■ ■ VqJU) , n = ^=4, with eigenvalues Eig(H) = - YT%=\ 2 cos (<?i + 2/-0- 
Let us define 

e ro =l-2 W ; ) m = 0,l,...,iV-l, (A.13) 

where the relation between g and m is that of (A. 9). We define m for m outside the interval 
in (A.13) by m = mm od7v- All the eigenstates of H indicated above are also eigenstates 
of O m with eigenvalues ±1, and it follows that 

N-l 

H = ©m cos (q(m) + 2/i) . (A.14) 

m=0 

The operators O m play a role analogous to the parameters €j for the eigenvalues of T(u): 
They fix the eigenvalues of H completely. Henceforth, we will consider them as constants, 
m G {+1,-1}, characterizing the eigenvalues of H rather than as operators acting on spin 
states. 



A.l Z 4 sector (N odd) 

By using the relation 

N 

JJcosx, = 2- (Ar - 1) sin (f + (-l)^>iV), (A. 15) 

i=i 

the expression for T(u, (i) of equation (A. 4) simplifies to 

N N 

T(u,fj) = ee -2 ^ expfi ^ e k x k ^ (cosw + % sinwe -2 ^ 3 J , (A. 16) 

fe=i j=i 

and the corresponding eigenvalue of Wd(O), found by setting w = 0, to 

v 

Eig(fi) =T(0,/i) = ee"— exp^^e^J. (A.17) 

k=l 

Because Ud(fi N ) = v Nd Ud(I), each eigenvalue of uid(fl) is v d times a f-independent phase. 
The dependence on v in (A.17) is hidden in the contribution exp(i Ylk=i e k x k)- By extracting 
from Xk the /^-dependent part and matching it with v d , one finds 

TV 

= d(-l)*P (A.18) 

fc=l 
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as announced in (3.16). This constraint reduces the number of possible eigenvalues to 
which is exactly the dimension of the representation Ud- 

To prove that the eigenvalues of Ud(T(u)) explore all the possibilities (i.e. that there 

are no repetitions) , we consider the XX Hamiltonian and initially show it can be written as 



H 



£((-iy + -e 



d+i ) sin 2x 



j- 



^odd, 



(A.19) 



$^((-iy + ^ 1 ©_iVi2 + 2 2 _i ( _ 1)J+ ^) sn,2.r y . 



N+d 



even. 



^ i=i 



This can be checked for all possible parities of 



%i and 



Here we only verify explicitly 



the case ^y 1 and both even. The other cases are very similar. Using the property 



m = 6 mmo dAf5 we thus compute 



N-l N-l 



H = J2 © m cos(^ + 2/x) = £ m sin(^±f^ + 2//) 

m=0 m=0 

iV-1 : ^T i JV-1 

= £e„» si .(iiffc + 2^(^ + £ )e„_^ S in(^ + 2^) 



n=0 

JV-1 
2 



E V + 2 ^)-E Bin(^ - 2„) 

n=0 n'=l 
2JV-1 2V-3 

E ©^sin(^ + 2/i)- ^ e_^sin(|^-2/i) 

r=l,5,9,... s=3,7,ll,... 
jV 

E(-l) J+1 e -f+^(-i)^ sin ( T + (-l) J+1 2^) (A.20) 



J=1 =2*, 



which is of the proposed form. 

Now, equation (A.19) can be directly related to the expression (A. 3) for the eigenvalues 
of UdiTi), provided the identification 



) ; ~ 1)J+ ^ -f + ^(-D-^' ^° dd ' 
e i ~ \ < iM+^n * 4 N+d (A.21) 

I {-iy+ 2 6 <m, ^even 



is made. As in Figure 4, the eigenvalues of i7 are non-degenerate (for generic /i), making 
this identification unique. This induces a bijection between eigenvalues of u)d(H) and of H. 



Because the [e=a) possible choices of m are all explored (exactly once) in the spectrum of 



36 




Figure 4: Z 4 sector (N odd). For N — 7, j — 1, . . . , N and = 0.02, the angles 2xj(/i = 0) 
and 2xj(fj), see (3.14), are indicated respectively in blue (full lines) and in red (dashed 
lines). Associated to a given angle labelled by ej, Q m is indicated in green (innermost) for 
even and orange (outermost) for odd. The groundstate for d — 1 is characterized by 
02, 63, O4 = —1 and by 03, ©4 = —1 for d = 3. In both cases, the remaining 6 m values are 
+1. For /1 generic, there are no degeneracies. 



H, the same is true for the Sj parameters associated with UdiH)- The proof of the selection 
rule for the overall sign e in (3.16) is presented in Appendix A. 3. 

In a given sector d, the identification of an eigenvalue as a groundstate energy depends 
on \i. To reproduce the selection rule given in [7], we consider // close to zero and study the 
largest eigenvalue of H. It is recalled that the eigenstates of H in the sector S z = ~ are of the 
form rj qi . . . V JU) with n = i.e. |{m; 6 m = -1}| = ^ and \{m; & m = = *±*. 
From equation (A. 14), the groundstate is built by choosing to be —1 the parameters 
6 m with cos(g(m)) negative and as large as possible in absolute value. It follows that the 
correct choice is 

f_i m = |M±2| ,...,N-1- 1^1, 

e m = \ ' L 4 J '. L 4 J (A.22) 

1+1, otherwise. 

(For fi outside the interval [—377, 277], this does not characterize the groundstate.) For ^f^- 
and ==i both even, for example, equation (A. 21) is rewritten as 

£2i = -Q aN+i _ v i = l,..., ^f 1 , e 2i+ i = { . N _, * N _, (A.23) 



4 



The negative m parameters are thus identified with the 6j parameters of the form e 2 i for i = 
. . . , It follows that e 2i+l = 1 for every i, and e 2i = 1 for % = . . . , ^f 1 -^. 

The only negative 6j parameters are e 2 i = e 2 i = — 1 for 2 = 1,..., here indicated in terms 
of barred and unbarred epsilons. This is precisely the statement of equation (3.17) for 
even. The other three cases (where ^tF^ or is odd) are examined in a similar fashion. 
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Comparison with the Z 4 sector selection rules of [7] : We now compare our selection 
rules for N odd with those of [7]. In that work, the selection rules pertain to the limit 
// — > and are expressed in terms of the excess parameters a and a defined in terms of our 
parameters €j and e,- = ejv+i-j by 



a 



a 



|{< = l,...,L^J;e* = -l}|-|{i = 0,... 1 L^J;e ai+1 = -l}|, 
|{i = l,...,L^lJ;e ai = -l}|-|{i = 0,...,L^J;c afrf i = -l}|. 



(A.24) 



Using (A. 18), we thus find 



a + a = \{i = e 2i = -1}\ - \{i = 0, . . . , ^1; e 2i+1 = -1}| 



2 ) '-M J \ I L ")•••) 2 ' 

JV-l V-l 
2 2 



iV 



I - E «*) - § - E = -t-i E(-!)^ 

i=l i=0 j=l 

d-1 <Z-1 



., • even. 

2 ' 2 



(A.25) 



which is how the corresponding selection rule is stated in [7]. It follows readily that a — a is 
an even integer. For the groundstate, from (A.24) and (3.17), we have 

f ^i, ^ even, 
o = e=\ ' * (A.26) 

^odd, 

again in accordance with [7]. 

A. 2 Ramond and Neveu-Schwarz sectors (AT even) 

First, we simplify the expressions for T(u,fi) by using the identities 

JL N /AT (sin/iiV, | odd, 

IJoosx^C-l^-^xJ (A.27) 
" [cos/ziv, | even. 

For both parities of |, we can thus write 

N N 

T{u,y) = ei% exp^y^efcXjfcj (cos w + i sin we" 2 "^ j . (A.28) 
fc=i j=i 

In the limit u — > 0, one finds the corresponding eigenvalue of u^(fi) to be given by 

N 

Eig(fi) = T(0, //) = exp (z J] e fc x fc ) (A.29) 

fe=i 
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and this must equal v d times a phase. The dependence on the winding parameter is 



exp(— i/jb €k), and hence ^2j=i e j = ~d as announced in (3.21). 



To show that there are no repetitions between the possible choices for the set of 



parameters {e,-} satisfying ^ 7 =i 



-d, we proceed as in the Z 4 sector by establishing 



a map between the eigenvalues of uj^T-L) and those of H. First, we note that the XX 
Hamiltonian can be written as 



H 



N 



3N_- sin 2xj, 



i=i 
N 

9 3N+2 ■ sin 2a; 



N 



even, 



3=1 
N 



3iv-2 j sin 2xj, 

3=1 



Y odd, | even, 



f odd, | odd. 



(A.30) 



For tt and | both even, for example, this follows from 



N-l 



N-l 



m=0 

N 



COS 



(2to+1)tt 
iV 



jv • sin 



(2i-l)7T 

N 



+ 2/i) = ^6 m sin((f -2m-l)f -2/i) 

m=0 

AT 

2^=-^e T _.sin(^-2^). 



(A.31) 



"V" 

= 2X, 



However, 



The other cases are examined in a similar way. 

Now, the final expression in (A.31) suggests the identification ej = — G>3al_, 
from the expression just before it, one might be tempted to make the alternative identification 
ej = This extra possibility is caused by degeneracies in the spectrum of u)d(H) that 

are not transferred to the spectrum of Ud(T(u)), see Figure 5. To resolve this issue, we recall 
that |{m; m = ±1}| = and note that with the identification €j = — 0w •, 



N 



N-l 



{m; T 



-1}| - |{m; 9 m = +1}| = -d, 



(A.32) 



m=0 



in accordance with (3.21). The same calculation based on the alternative identification, 
€j = yields ^ . e j = d and is therefore incompatible with (3.21) for d > 0. As 

discussed in Appendix A. 3, even for d = 0, only the first identification results in a fully 
consistent picture in general. 

As in Appendix A.l, the identification is a one-to-one map between the eigenvalues 
of the two models. Because all the different choices of Q m are explored in the XX model 
(with the constraint |{m; G m = — 1}| = —5-^), there are likewise no repetitions in the column 
configurations characterizing the eigenvalues in the loop model. In general, the identification 
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-e 3 -e 4 



-e 4 




(a) (b) 

Figure 5: For N = 8, j = 1, . . . , N and \i = 0.04, the angles 2xj(fi = 0) and 2xj(fi), see (3.14), 
are indicated respectively in blue (full lines) and in red (dashed lines), in (a) the Ramond 
sector (N even, | even) and (b) the Neveu-Schwarz sector (N even, | odd). The groundstate 
for d = is characterized by 2 , ©3, ©4, ©5 = —1 and the one for d = 2 by 3 , 4 , 5 = — 1, 
with all other 6 m values equal to +1. In both sectors, because sin2a; ? - = — sin2x,-, n, the 
eigenvalues of H are degenerate and, a priori, this could lead to two possible identifications. 

depends on the parities of y and |, as we have 

N 

j even, 

f odd, I even, (A.33) 
f odd, I odd, 

in accordance with (A. 30). 

The proof of the selection rule for the overall sign e in (3.21) is presented in Appendix A. 3. 
As in Appendix A.l, the groundstate of H is obtained by fixing the parameters m as in 
(A. 22). Here this translates into the selection rule (3.22) when (A.33) is applied. 

Comparison with the Ramond sector selection rules of [7]: We now compare our 
selection rules for iV and | both even with those of [7]. In that work, the selection rules 
pertain to the limit fi — > and are expressed in terms of the excess parameters a and a 
defined in terms of our parameters 6j by 

o = |{j = l,...,[^J;e, = -l}|-|{ J = l,...,L^J;e^ +1 _, = +l}|, 

(A.34) 

a = \{j = i, . . ., LfJ; e, +1 _ 3 = -i}| - |o = 1 LtJ; e f +j = 



-Q 3N 

4 



-3 ' 



-0 3JV-2 ., 
4 J 
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Using J2f=i e j = we thus find 

- + * = \ ( m - L E - \ ( m + L E 

LfJ LfJ 

+ !(LfJ-E e f + i-i)-l(LfJ + E e 



L^J \S+2l\ 



3=1 3=1 
N 

= -iE* 

= 5 (A.35) 

which is the corresponding selection rule of [7]. It follows readily that a — a is an even 
integer. For the groundstate, from (A. 34) and (3.22), we have 

o = ° = ~, (A.36) 

again in accordance with [7]. 

Comparison with the Neveu-Schwarz sector selection rules of [7] : We now compare 
our selection rules for N even and | odd with those of [7]. In that work, the selection rules 
pertain to the limit /i — > and are expressed in terms of the excess parameters a and a 
defined in terms of our parameters tj by 

a = -8 + \{j = 1, . . . , [f J; e 5 = -1}| - \{j = 1, . . . , [f J ; e N -j = +1}|, 

(A.37) 

* = 10' = i,-,L i ¥J;^- i = -i}|- 10' = i,---, L J ¥J; e f + i = + 1 >l- 

Here we have introduced the two additional parameters 

8,8 €{0,1} (A.38) 

to implement the freedom in the definition in [7] of the excess parameters. It is also noted 
that €n_ and do not appear in the expressions for a and a. We now evaluate 

2 

LfJ LfJ 



= -*-*+l(L?J-X>H(LIJ +!>-/) 

3=1 3=1 

L^J L^J 

+ l(L J ¥J-E e f-i)-l(L i ¥J+E e 



2 

i=i i=i 

JV 



JV +i 



^ + 1(2 + 6^ + 6^-25-25) (A.39) 
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where we have used Ylj=i e j = ~d. We thus recover the selection rule 



J 3 



of [7] by setting 



d-2 , 
a + a = — — (A.40) 



5=^—^, 5 = S (A.41) 



in which case 

a = = LfJ; e 3 = -1}\ - \{j = 0, . . . , LfJ ; e N - S = +1}|, 

(A.42) 

* = Kj = i,---, L^J; = -i>| - |{j = o,..., L^J; = +i}|. 

It follows readily from the selection rule (A.40) that o — a is an even integer. For the 
groundstate, from (A. 43) and (3.22), we have 

a = a = — , (A.43) 

again in accordance with [7]. 
A. 3 The overall sign e 

We now prove the selection rules for e given in equations (3.16) and (3.21), again by resorting 
to the connection with the XX model. As shown in [12], the representative of in the 
XX representation of £PTL]y(a, = 0) is given by the matrix 

fi±i = r,±^"t±i. (A.44) 

The operators t and t _1 introduced here translate the N spins to the left and right by one 
position, implying that 

t^ap* 1 = a a m . (A.45) 

For both parities of N, every eigenvalue of ujd(T{u)) (and u)d(H)) is associated to an 
eigenvalue of Ud(Q), see equations (A. 17) and (A. 29). It is shown in [12] that every element 
of SPTLn{c(, P = 0) has the same spectrum in the uid representation as in the XX represen- 
tation with S z = |. This is therefore true for T(u), H and f^ 1 , in particular. To determine 
e, we use the identifications (A. 21) and (A. 33) along with the equality Eig(fi) = Eig(0). 

N—d 



First, we seek the eigenvalue of f2 associated to the state rj qi . . . rj qn \U), where 



n 



2 • 

The operator v 2S * acts on this state as the multiplicative constant v d . By applying the 
translation operators to the fermionic operators in Proposition A.l, one finds 

r x Cjt = (-af)cj+i = c j+1 (-al), j = 1, . . . , N - 1, 

r x c N t = (-^) Ci (-1)t = Cl (_^)(_i)f +S' (A . 46) 



42 



and 

r\t = (-oJ)(e% + *=( e -* N (-l)% +s ' +1 - 1)) (A.47) 

= (/% + -^( e -^(-l)f +^ - (A.48) 

Applying t^ 1 on eigenvectors of H yields 

t-\, ■ --vju) = (rV) • • • (rV)*' 1 ^ = (-i)v^*^ . . .^jc/> (A.49) 

where we have used the two expressions (A.47) and (A.48) to rewrite t~ lr q q .t for j even and 
odd, respectively. Because 

J m le m =-i m|e m =-i He m =i ™=0 

1^ 7T f AT — 1 ^odd, 

= --Ve m g(m) + -x ^ (A.50) 
2 ^ mHK J 2 ^ ^+d e ven 

m=0 I > 2 cvcli ) 

one readily finds 

• N ^ 1 f i N + d oc [d 

Eig(fi) = ^(-O^expQ 8 m ?(m)), a = ' ^ ' (A.51) 



?Ti=0 



The next and final step is based on the identifications (A. 21) and (A. 33) of the 6j and 
parameters, and is performed separately for the two parities of N. 

1*4 sector (N odd): The selection rule for e we wish to establish states that 

— 1) * , V 1 even, 

d J ( A -52) 

i \ d + 1 rf— 1 l l v ' 

-1) 4 , odd. 

Let us be specific and consider — ^ and both even. In this case, the identification (A. 21) 
specializes to (A. 23) which is equivalent to 



N-l 

T ^+2m) if* — u, ... , 4 , 



ejv+i, 9rri , m = 0, . . . 



6 m = <( -e3^i_ 2m , m = ^,...,^, (A.53) 
e 2m _3^, m=^tl,...,JV-l. 
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To evaluate Eig(Q) in (A. 51), we now compute 

/.N-l \ A V ^ JV-l 

exp ^^O m g(m) = exp ^ ( J] meiv±i +2m - J] me3^ti_ 2m + ^ me 2m „^_i 

V m=0 / \ m=0 ' m= N + 3 m= 3iV+i 

iV-l N-l N-5 

= ex P ( E - ¥fe + E " + E (2^3 + ^Ka+i)] 

V ji=^fi J2=l J3=0 

, V-l -V-l V-5 

^ E^' - ^ J ex P [j ( - E £ 5Wi+i - 3 E t2 n + 3 E £ 2J3+i) 



JV-l 



AT \ / „. JV 



ex p bjvE( 2j ' " ^3 ex p _ jD -1 )^ k p E 

V i=i J \ k=i / \ ii=^i 



0=1 

where we have used (A. 18). This allows us to compare Eig(f2) with 

N 

ESg(n) = ee-^^exp(^^(2j - l) CjF ), (A.55) 

3=1 

from which we deduce that e = (— l) - * 1 . Similar arguments carry through for the other 
parities of and 2=A ( thereby completing the proof of the selection rule (A. 52). 

Ramond and Neveu-Schwarz sectors (iV even): The selection rule for e we wish to 
show is given by 



-1)4 7j even, 

iV^ d j j (A.56) 
-1) 4 , I odd. 



e 



For f and f both even, for example, the identification (A. 33) is equivalently given by 

-63N _ m , m = 0,...,^-l, 



Q m={ / 3N \ T , (A.57) 



-em_ m , m = ^,...,iV-l, 
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and one computes 



/ . N-1 \ / . "4~ 1 N-1 

eXP U E @m ^ = eXP W ( ~ E ( 2 ™ + 1 ) £ ^-m ~ E + ^^-m) 

\ m=0 / \ m=0 m = ML / 



3N 

— N 



i=l i= ML+l 
TV TV TV 



— E(2j-ih-— E e i- m E e ^ 

= «p(^ffi( - D%) exp(^)(-l)f (A.58) 
By comparing Eig(f2) with 

TV 

Eig(tt) = ez Vex P (^E(^' " ^ (A.59) 

3=1 

valid for y and | both even, we find e = (— 1)*. Similar arguments show that the selection 
rule (A. 56) also holds for the other parities of ~ and |. 

The calculations above and the ensuing comparison with Eig(fi) could have been carried 
out for the alternative identification 6j = On - discussed following (A. 31). In that case, we 
would have obtained 

TV 

6 = (-1)* exp(- ^ E( 2 ^' - V**) > ( A - 6 °) 

but this is in general incompatible with the requirement e£ {+1,-1}. Thus, even in the 
d = sector, only the first identification yields a fully consistent picture in general. Similar 
arguments apply for the other parities of y and |, thus demonstrating that the identification 
(A. 33) is the correct one in the general case. 



B Loop model on helical tori 

Here we generalize the loop model discussed in the bulk part of this paper by considering the 
model defined on helical tori [31-33]. Prior to forming the torus, the loop model is defined 
on a finite lattice, which we represent diagrammatically by a rectangular array of M x jV 
square tiles. In this planar representation, the toroidal boundary conditions are encoded 
in the two periodicity vectors i/i, i/ 2 £ As indicated in Figure 6, we parameterize these 
vectors as 

z/ 1 = (iV,0), z/ 2 = (t,M), (B.l) 
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where t is the helicity. Its role is as follows. To construct the torus, one first identifies 
the right and left edges of the rectangle such that the tiles cover the outer surface of the 
resulting vertical cylinder. A torus is then formed by gluing together the upper and lower 
edges of the cylinder, but before doing that, one twists the upper edge of the cylinder t tiles 
in the clockwise direction seen from above. The geometry of the helical lattice torus is thus 
defined by the triple (M,N,t), where t is negative if the twist is in the counterclockwise 
direction. In the bulk of this paper, we have exclusively considered t — 0, but the results 
can be generalized in a straightforward manner as we will discuss in the following. 




Figure 6: (a) A configuration c of the loop model on a 4 x 6 torus with helicity t — 2, with 
five contractible loops and two non-contractible loops of homotopy {0, 1}. (b) The weight 
J-2(c) of the configuration c is equal to T{c ■ Q~ 2 ). 

Loops (and similarly clusters) of homotopy {a, b} now wind around the torus a times in 
the vi direction and b times in the v% direction, with the same convention for the sign of a 
as in Section 2.1. Cross-topology clusters are those that wrap the torus along both v\ and 
Vi- Of course, in the FK model, we restrict the helicity t to even values only. 

For a fixed connectivity c, the homotopy of non-contractible loops and the number of loops 
in general depend on t. The action of the abstract functional J-(c) that assigns the correct 
weights to configurations when the cylinder is closed into a torus (defined in equation (2.37) 
for t — 0) depends on the helicity and is here denoted by •T-i(c) (with J-q = J-). However, 
as illustrated in Figure 6 (b), studying a torus with helicity t is equivalent to considering 
connectivities times Q~ l on a torus with trivial helicity: 



To calculate the partition function Z£, our prescription (2.42) for T{c) needs to be slightly 
modified: Connectivities contributing to c = T M (u)il~ l can have up to M + t loops crossing 
the virtual boundary. This suggests changing the upper and lower bounds of the sum in 



(B.2) 



g d (fi,a) to ±(M + t) 



M+t 




(d > 0). 



(B.3) 



k=-(M+t) 
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This simple prescription is not optimal, though, since it may include vanishing terms. Indeed, 
because the first t layers of boxes are just some terms, the one corresponding to k = M+t 
for instance, are simply zero. Nevertheless, with this slight modification of the definition of 
Qd(fJ>,oc), equation (2.42) remains valid. 

The problem of computing the partition function Z\ therefore boils down to calculating 
the eigenvalues of Ud{T M {u)Vr l \ The eigenvalues of u^(fi) corresponding to T(u,fJ,) are 
given in equations (A. 17) and (A. 29), so in general 



Eig(T M ( M )fi-') = (KQm)) 



iirNt 



M M—t 

x e x e 4 



(_l)iV+l 



n 

3=1 



g ite j x j 



(e m + 



ie iu ejtanxj) M (B.4) 



and the generating functions Gd{z, v) should be modified in a similar way to depend on the 
helicity t. Below, we write G d (z,v) for these modified generating functions and G d (z,v) for 
their normalized versions defined as in (3.40). 

Finitizations of the generating functions G d (z, v) are obtained by looking at finite exci- 
tations for M,N,t 3> 1, not only for fixed aspect ratio (3.42), but also while keeping the 
ratio 

7 = ± (B.5) 
fixed. With the parameterization (3.41) for /i, elementary excitations are given by 



lim 

M,JV,t>l 
M=8N, t=jN 



e m - ie~ lu - 



e iu _|_ %u £ an x 



tanxA M _ E{a) 



(B.6) 



where 



lim e 

M,JV,t>l 
M=SN, t=-yN 



2il..r 



fU) 



e — ie 



+ ie 



tanx m \M _ 
tanx m ) 



-l) M+t q^~ a \ 



(B.7) 



EM) 



and 



2.7-1 
4 



2j-l 
2 



+ ("!) 



3 + 



d+l 

2 a, 



N odd, 

A" even, | odd, 

d . 

*> 2 



N+l-j 
f + 1-3, 



N odd, 
A" even 
N even 



' 2 
d 
' 2 



odd, 
even, 



g 7 = exp ( - 27ri(5e- 2iu - 7)) , g 7 = exp (2vrz(5e 2m - 7)) . 



(B.8) 
(B.9) 



The expressions for G a (z, v) are now obtained from (3.50), (3.56) and (3.57) by replacing q, q 
by 9 7 , q-y and (— 1) M by (— l) M+t . For A" even, because the eigenvalues of Wd(fi) corresponding 
to the groundstates in the d = and d = 2 sectors are equal to 1, the ratio (3.58) of the 
maximal eigenvalues of ^(T^fi - *) in these sectors is still given by (qxqx)^ = (qq)* in the 
continuum scaling limit. 

In this limit, for M even, the partition function for critical dense polymers on the torus 
with helicity t = 'jN is then given by (3.66) or (3.67) for A" odd and even, respectively, 
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with the 7-dependent definition (B.9) for the nome. For M, N and t all even, this partition 
function is therefore modular invariant. 

As discussed in [38] on the related lattice spin models, the spectral parameter (for < 
u < |) can be interpreted as measuring the spatial anisotropy of the lattice through the 
anisotropy angle 

7T?/ 

6 = = 2u. (B.10) 
A 

With this geometric interpretation, changing the spectral parameter away from the isotropic 
point u = -| = j thus corresponds to distorting the elementary square faces (1.1) on the 
lattice into rhombi whose bottom-right angle is given by 9. It follows from the discussion 
above that a nontrivial helicity has the effect of changing the spatial anisotropy of the lattice, 
in accordance with the situation for finite lattices as in Figure 6 (a). Concretely, the helicity 
acts on the modular parameter (see (3.45)) 

T = 5e i{ *- e) = -5e~ 2iu (B.ll) 

as 

T -). T ' = r + 7. (B.12) 
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